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j> ! Abstract 

O; 

The proton's transversity distribution will be measured at BNL's Relativistic Heavy Ion Collider 



in upcoming experiments using the transverse Drell-Yan process. Understanding the one-loop cor- 
rections is therefore important. Here, the collinear structure in transverse Drell-Yan is investigated 
in detail using dimensional regularization and the correct behaviour is found, although the mech- 

^ anism is non-trivial. The resulting n-dimensional transversity splitting function (and consequently 

the one-loop transversity distribution and its two-loop evolution) is found to be the same in both 
the anticommuting-75 scheme and the HVBM scheme. Alternative schemes are considered. 
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1 Introduction 



Recently, two-loop calculations [1,2,3] of the transversity splitting function A T P qq [4] relevant 
to the evolution of the transversity distribution A T f q (or h\) [5] (to be measured in upcoming 
experiments at BNL's Relativistic Heavy Ion Collider) have been carried out using dimen- 
sional regularization (DREG). Of some concern is the fact that an anticommuting 75 was used 
in the above determinations. This is known to be mathematically inconsistent. Fortunately 
[2] the traces only involve even numbers of 75's, so we do not anticipate any inconsistencies 
for this specific determination. Of greater concern is whether DREG itself is suitable for the 
calculation of higher order corrections to processes with transverse polarization, either in the 
anticommuting-75 scheme or the mathematically consistent 't Hooft-Veltman-Breitenlohner- 
Maison (HVBM) [6,7] scheme. More precisely, it is necessary to verify the correct behaviour 
of the squared amplitude, relevant to some subprocess involving transverse polarization, un- 
der collinear gluon bremsstrahlung. That behaviour must be consistent with that of other 
regularization schemes in order to yield a meaningful and process independent Axfq- This 
is the issue we shall address in this paper. 

2 Generalities 

The process we shall consider is transverse Drell-Yan, where we have 

A(P 1} Si) + B(P 2 , ±s 2 ) -> r( P3 ) + l + (p 4 ) + X. (1) 

Here A, B denote hadrons with momenta Pi, P 2 and transverse spin vectors Si, S 2 . The spin 
vectors satisfy Si ■ Pj — 0, Sf — — 1. In leading order, this process is mediated by 

q(Pi, si) + qfa, ±s 2 ) -> 7 *(g) -> Z"(p 3 ) + 1 + (Pa) + (q <-> q), (2) 

shown in Fig. 1. In the parton model pi = X\Pi, p 2 = x 2 P 2 for some < Xi,x 2 < 1. At 
next-to-leading order in QCD, we have the gluon loop corrections to the above subprocess 
and the bremsstrahlung subprocess(es) 

q(Pi, si) + q(p2, ±s 2 ) -> 7 *(g) + g(k) -> Z"(p 3 ) + l + (p*) + g(k) + (q <-> g), (3) 

shown in Fig. 2. There is no qg subprocess contributing to transverse Drell-Yan. For the 
extraction of Axfq, one is interested in the transversely polarized cross section and the 
corresponding subprocess cross section defined by 

A T a = l -[a{Si, S 2 ) - a(Si, -S 2 )\, A T a = ±[a(si, s 2 ) - a(s u -s 2 )\, (4) 

respectively, with Si — Si. We see then that the Si (sj) represent the "up" directions, or spin 
quantization axes. A T a is obtained by appropriately convoluting the A T f q (xi) with A T a. 
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Comparison with experiment then yields information on Ay/ ? . We may define the invariants 
M 2 = q 2 = sw, s = (pi + P2) 2 , t=(p 1 - p 3 ) 2 , u=(p 2 - p 3 ) 2 . (5) 



If we let the component of p 3 transverse to the beam axis (z) define the x axis, then in 
n = 4 — 2e dimensions, in the pi,p2 cm., the momenta and spin vectors are decomposed as 

y x z 

si,2 = (0;0,...,0, sin 0i j2 , cos 0i i2 , ), 

p 3 = |p 3 |(l; 0, . . . , 0, sin 3 , cos9 3 ), Pi,2 = ^(1; 0, . . . , 0, ±1), 

fc = |fc|(l; . . . , sin 9 ljk sin 9 2 , k cos 9 3ik , sm 9 hk cos 9 2 , k , cos 9 ljk ), (6) 
where, for fixed M 2 and ^3 direction (but not magnitude), 



\P 



M 
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2 [% + I^Ksin^sin^i^cos^jfc + cos6> 3 cos 6 hk )} ' 



\k\ = ^(l-w), q = y^(l + w). (7) 
In k, the . . . represent the n — 4 components to be (trivially) integrated over. 

To obtain a nonvanishing result, we must not integrate fully over the azimuthal angle of p 3 , 
3 . Writing d n ~ l p 3 = \p 3 \ n ~ 2 d\p 3 \d n ~ 2 Q 3 , we may present the relevant 2^3 phase space, 



A T da 1 M 2 ~ i£ I A; I 1 " 26 



dM 2 d n ~ 2 Q 3 (2tt) 5 - 4£ 2 2 ~ 2e s 3 / 2 

/< 



(8) 



x n — = ™ . 

[q + I A; I (sin 6 3 sin 6*1^ cos #2^ + cos 9 3 cos #i,fc)] 2 ~ 2£ 



where \M\ 2 is the color averaged 2 — > 3 particle squared amplitude. In general [8] 



2tt 



Jd m ' 2 n = J d9 1 sin m - 3 9 X J d6 2 sin m ~ 4 9 2 ...j d9 m _ 3 sin m _ 3 J d9 m _ 2 . (9) 



The 2 — * 2 particle phase space is simply 

A T da 2 2£ M~ A - 2£ 
dM 2 d n ~ 2 Q 3 ~ (2tt) 2 - 2£ 



5(1 - w)A T |M| 2 _ >2 . (10) 



Note, our squared amplitude normalization corresponds to the convention J2x U (.P> ^) u (p, A) = 
2V 2 for p 2 = 0, with A denoting helicity. Once the above differential cross section contribu- 
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tions are obtained and the virtual contributions and factorization counterterms are added, 
yielding a finite result in the limit e — > 0, one can integrate over # 3 and obtain A^do j dM 2 d^. 
The original calculation of the latter quantity, at one-loop, was done in [9] using the regu- 
larization scheme where the gluon is given a mass in order to control the collinear and soft 
singularities. That scheme had already been used successfully in unpolarized Drell-Yan [10]. 
Then, in [11] the calculation was done using regularization by dimensional reduction (DRED) 
[12]. The general approach for converting from DRED to DREG, for arbitrary polarization, 
was given in [13]. That approach requires knowledge of the n-dimensional transversity split- 
ting function, and a more careful derivation of it was necessary. This was done in [2] , where 
the results of [9] were converted to DREG. Unfortunately, the n-dimensional transversity 
splitting function was explicitly derived in [2] only in the context of a two-loop calculation. 
No details regarding the collinear limit were given, other than that the correct form resulted. 
The result of [2] for the transverse Drell-Yan cross section confirms the earlier DRED result of 
[11] and it confirms the general form of the cross section given in [13], valid for all consistent 
n-dimensional regularization schemes - the scheme dependence is correctly parametrized in 
terms of the n-dimensional splitting function. 



3 The squared amplitude 



Let us write for the 2-^3 particle amplitude M 
Figures 2(a),(b) respectively. Then 



Mi + M 2 , where M 1 , M 2 are shown in 



A T |M| 2 = A T |Mi| 2 + A T |M 2 | 2 + 2A,-.U I .U,' 
1 



-C< 



+ 



Tr[7s fa Al^i- fa Alvitfi- 



[(pi - k) 2 } 2 

Tr[ 75 fa i><fi^ 2 - #)7a75 fa Alpfa- 



1(P2 - k) 2 ] 2 
2 



(pi - k) 2 (p 2 - k) 

( M 2 N 



r Tr[ 75 fa Alai^i- /07^75 fa Alpfa- 



where C is an overall factor. Let Ar|M^| 2 denote the n-dimensional Born term in DREG. 
Writing At|M^| 2 = At\M^\ 2 + sAt\M%\ 2 , we have in the anticommuting-75 scheme, 



A T \M n B \i c = C B 
= C B 



2M p 3 ■ si p 3 ■ s 2 + tusi ■ s 2 + £—— si ■ s 2 



M 



M 4 



! sin 6» 3 cos(0i + <p 2 - 20 3 ) - e— cos(0i - 2 ) 



(12) 
(13) 
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where Cb is an overall factor and clearly w — 1. We have now rotated back to a frame with 
an arbitrary x axis direction via lj2 — ■> 0i,2 = 0i,2 ~ 03- In the HVBM scheme 

Ar|M||| VBM = -A r |M||l c . (14) 

At this point, we notice that the Born term has the wrong azimuthal dependence in DREG. 
Since the effect is order e, it will not manifest itself as long as the Born term factors properly. 
Then the e-dependent part will cancel along with the singularities that multiply it. This is 
the case for the virtual corrections and the soft corrections pose no problems thanks to the 
Bloch-Nordsieck mechanism. The remaining question, therefore, is: what is the behaviour in 
the limit of collinear gluon bremsstrahlung? 

4 The collinear limit 

In order to investigate the collinear limit, k || pi, we define the quantities 

A T P< q 4 (z) = C F (j?- - 2) , A T P qq (z) = -C F (1 - z), (15) 

where the superscript < indicates that z < 1. Let us also define 

t' = (p'i-p' 3 ) 2 , W = (p' 2 -p' 3 ) 2 , p[=w Pl , p' 2 =p 2 . (16) 



p' 3 is the value which p 3 takes in the above collinear limit, for fixed Q3 and M 2 . Hence 
the primed quantities are those relevant for the kinematics of the Born term which should 
factor out in the above limit. Now, in the collinear limit, k \\ p±, the 2 — > 3 particle squared 
amplitude in the anticommuting-75 scheme takes the form (with w < 1) 



At|m| c 2 o1 , ac 



4C 



wpi ■ k 

— (ewtu + 



AT|MgUc(i - f -W<.v )+ ,A T p M w] 



CbCf 
e 2 M 4 k- Sl k- s 2 
2 ' p x ■ k 



,, 2 k ■ sAwt + u + M 2 ) 

+ eM 2 p 3 ■ s 2 ; 

Pi ■ k 



(17) 



plus terms which are nonsingular and do not give rise to scheme dependences - the 0(s 2 ) 
terms above may also be dropped. In the HVBM scheme, we find 



At|M| c 2 o11jHVBM = A T |M| c 2 olljAC (e <- -e) + 0(e 2 , ek 2 ), 



(18) 



so that the structure is identical in both schemes, with only the sign of e reversed (including 
in the Born term). There are no finite contributions from k integrations, where k is the 
vector whose only nonzero components are the components of k having index greater than 
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three (the timelike index being zero) in n dimensions; they contribute like an extra power of 
e. Also, none of the 0(e 2 ,ek) terms dropped multiply a soft divergent term. 

The term ~ k ■ s 1 k ■ s 2 in (17) depends on the "azimuthal" angles # 2 ,fc and 6> 3ifc , of the gluon, 
in n dimensions. In the collinear limit, k \\ p\ (i.e. 9i ik —> 0) those angles must be integrated 
over since they are unconstrained. From (7), (8), we see that the azimuthal dependent part 
in the denominator of the phase space and that of t, u vanishes in the collinear limit by order 
sin^fc (see below), so that the only finite dependence of that term on 02,k,&3,k comes from 
the factor ~ k • s± k ■ s 2 . This term gives a finite contribution due to the 1/pi ■ k factor. Since 
it is multiplied by e, it only gives a finite contribution to the cross section from the phase 
space region 9i jk — > 0, where the 1/e pole coming from integrating the overall 1/pi ■ k factor 
arises. Hence, we may perform the azimuthal integration in that limit. From (6), (8), (9) we 
see that we may make the effective substitution 



k ■ s\k • s 2 
Pi ■ k 2 



Pi ■ k 



Pi ■ k 



TV 71 

J de 2 ^ 2e e 2)k J de 3ykS m^- 2£ e^ k 1 



1 — w 



cos(0i - (j) 2 ) 
(1 - e)p l ■ k 



w) 



si ■ s 2 



1 — e p\ ■ k 



(19) 



Note, it would make no sense to retain the 0(e) term since our approach only determines the 
finite contribution. Similarly, the last term of (17) picks up an additional finite azimuthal 
dependence from the azimuthal dependence of t and u, which is of order sin 9i tk and can be 
obtained by series expanding (7) about sin^x ^. = 0. Without this extra dependence, the term 
would not contribute. The procedure is then similar to that used for the term ~ k ■ s 1 k ■ s 2 . 
We end up with the effective substitution 



k ■ Si(wt + u + M 2 ) 2(1 -w)s 1 -p 3 



Pi ■ k'< 



Pi ■ k 



+ 0{e). 



(20) 



Substituting (19) and (20) in (17), we see the explicit cancellation of the term ~ A T P qq (w): 

AC A T \M%\ 2 (t = t') 



At|M| c 2 o11 = 

wpi ■ k 



CbCf 



A T P<Aw) + 0(e 2 ), 



(21) 



for both the anticommuting-75 scheme and for the HVBM scheme. This demonstrates the 
required factorization property in the collinear limit. We also confirm the finding of [2] that, 

A T P< q n (z) = A T P< q \z), z < 1, (22) 
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in dimensional regularization for the anticommuting-75 scheme. In addition, we have shown 
that ATP< ,n (z) is the same in the HVBM scheme. The 5-function part has an e-dependence 
equal to that of the unpolarized P q {z). According to the "+" prescription [14] one obtains 



As this result agrees with the one obtained as an intermediate step in the determination 
of the two-loop AxP qq [2], in the anticommuting-75 scheme, which was not specific to any 
particular process, this provides a good check of the process independence of A T P™ q {z), and 
hence of A T f q . 

In order for A^/^dreg to be meaningful, we should be able to relate it to A T f q in some 
other scheme. It is easy to see that terms ~ l/pi ■ k 2 in (17) cannot arise in four- dimensional 
schemes with massless quarks - by four-dimensional, we mean four-dimensional tensor struc- 
ture. Terms ~ 1/pi-k 2 , ~ l/p 2 -k 2 can only come from A T |Mi| 2 , A T \M 2 \ 2 respectively. These 
squared amplitudes vanish in a trivial fashion in four dimensions, so that only AtMiM^ 
contributes. In n-dimensions this is not true. It is violated by order e, due to the relation 
Itilplal^ — ^9pal — 2e7 p 7 CT . The most natural four-dimensional alternative to DREG is 
DRED, where the tensors and gamma-matrices are kept four-dimensional, but the momenta 
are taken to be n-dimensional with, formally, n < 4. In that scheme one obtains the collinear 
behaviour as in (17), but with e — > everywhere, including in the Born term factor. Hence 
the collinear structure is the same as in (21), except with the four- dimensional Born term. 
Since DRED and DREG have the correct factorized form in the collinear limit, we can use 
the technique of [13] to convert subprocess cross sections from one scheme to another. Differ- 
ences arising in the subprocess cross sections are cancelled by differences in the transversity 
distributions, calculable using the various A T P qq . The conversion formula is given in [13] and 
it has the same form as in the unpolarized and longitudinally polarized cases (see also [15]). 

It was checked in [15] that the differences in the two-loop evolutions of the longitudinally 
polarized parton distributions could be traced back to the differences in the n-dimensional 
one-loop longitudinally polarized splitting functions (as well as the differences in one-loop 
factorization scheme), or equivalently, the differences in the one-loop parton distributions 
themselves. Since there are no such differences between the A^-P™ of the anticommuting-75 
scheme and that of the HVBM scheme (meaning the transversity distributions are also the 
same at one-loop), we conclude that the two-loop evolution of Arf q , and consequently the 
corresponding two- loop A T P qq , should be the same in both schemes as may be checked. 

5 Alternative regularization schemes 

We now consider the details of certain alternatives to DREG. Although not specifically 
related to transversity, this is a good opportunity to investigate alternatives such as DRED 
since DRED is free of the many complications we had to deal with in DREG in the study 
of transverse Drell-Yan. Hence its usefulness may be appreciated here. As was pointed out 
in the original one-loop DRED calculation of the transverse Drell-Yan process [11], a UV 



A T PVz) = C F 



2 



2 + d + l)5(l-z) 



(23) 
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counterterm must be added to the vertex correction (see also [13]). This is not a problem 
since the counterterm may be generated by a process independent Feynman rule. Similar 
terms were pointed out in [16]. Still, it is useful to consider other alternatives. 

For Drell-Yan (or deep-inelastic scattering) the solution is simple. One simply uses the n- 
dimensional metric tensor in the virtual photon propagator. Consequently, the g a @ in (11) 
becomes the n-dimensional one, whereas the g^ u and the gamma matrices remain four- 
dimensional. This projects out only the physical part of the vertex loop. Then, one also 
finds the correct behaviour in all the collinear limits, including for the unpolarized and lon- 
gitudinally polarized cases and for the qg subprocess which also arises there - one finds the 
relevant four- dimensional splitting function multiplied by the n-dimensional Born term. The 
n-dimensional Born term which factors out is the same as that of the anticommuting-75 
scheme of DREG. This is not surprising if we note that by taking all metric tensors aris- 
ing from bosons (including qg vertices, etc. . . ) to be n(< 4)-dimensional, DRED reduces to 
DREG (in the absence of the 75 problem), up to a possible DRED mathematical inconsis- 
tency, which only rarely arises. Hence the Born terms are the same. The approach of taking 
all metric tensors arising from bosons to be n-dimensional defines a scheme, which we shall 
denote as scheme II (scheme I being DRED with counterterms) . In scheme II, we observe a 
tradeoff between the ill-definedness of an anticommuting 75 (and analogous e^^) in n > 4 
dimensions and that of g^ u in n < 4 dimensions. Whatever the approach used, in n < 4 
dimensions one must first contract all inately four- dimensional tensors such as e^ vpu with 
each other before contracting with the n-dimensional metric tensor as the contraction of the 
above quantities is ill-defined [17]. For situations involving traces with an odd number of 
75's, or explicit occurrences of the £-tensor, scheme II suffers from the same arbitrariness 
as the anticommuting-75 scheme of DREG. The advantage of anticommuting-75 schemes 
like scheme II is that they generally satisfy Ward identities relevant to electroweak interac- 
tions (and thus require no counterterms) whereas the HVBM scheme generally requires UV 
counterterms to restore those identities as well as various finite renormalizations. 

We now return to the first variant on DRED discussed above; the statements below are not 
relevant to scheme II. It is permissible to take the virtual photon as being n-dimensional, but 
the radiated gluon must remain four-dimensional for consistency (i.e. same splitting functions 
and UV sector) with the usual DRED. As well, the latter directly leads to the vanishing of 
A T |Mi| 2 , A T \M 2 \ 2 as can be seen from (11). The above approach may be stated as a more 
general rule for DRED which so far appears to be valid at one-loopQ Metric tensors arising 
in propagators of bosons not occurring in virtual loops and which ( unless massive ) never go 
on shell (where a massless boson might develop a mass singularity or a soft divergence) are 
taken as n-dimensional. This saves us from having to add counterterms in many cases, but 
not in all cases. Initial and final state boson lines are taken as four- dimensional in this scheme. 
Formally, one should add counterterms in DRED, so in our case the above procedure simply 
amounts to a trick to avoid that. One can always add the counterterms in this scheme. Then 
many of them will simply decouple. If conventional DRED should fail in some circumstance, 
it may prove useful to make the above approach, which we denote as scheme III, more formal, 



2 We have not explicitly checked what happens when 3-gluon vertices are present. 
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since it could fix problems which counterterms are not capable of fixing. One such case which 
motivates further study of scheme III (and also II) is the seeming problem with conventional 
DRED pointed out in [18]. Since that problem is not in the UV or soft sectors, one does not 
expect the usual counterterm approach to work. 

For the longitudinally polarized qg subprocess, in scheme III, all four-dimensional algebra 
including the e-tensor contraction was performed first, then the contraction with g a/3 was 
carried out. This simple approach leads to incorrect results in scheme II. Additional ad-hoc 
rules must be introduced. One prescription considered was to first take all the traces in four 
dimensions (by using n-dimensional metric tensors to project indices occurring in the traces) 
and contract out all resulting repeated indices, then replace all remaining n-dimensional 
metric tensors with four dimensional ones, contract out the e-tensors (so far untouched) and 
perform all the remaining contractions in four-dimensions. Then, the correct form results, 
with the four-dimensional splitting function times the n-dimensional Born term arising in the 
collinear limit. The prescription may be stated as follows: The e-tensors must be contracted 
at the end, with only four- dimensional metric tensors remaining. Those metric tensors are 
obtained by replacing all remaining n-dimensional metric tensors with four dimensional ones 
(i.e. gP a — > ) after all repeated indices have been contracted out. It should be noted that 
for traces involving an even number of 75's, the possible mathematical inconsistency is easily 
removed by performing such traces in n dimensions. Whether n < 4 or n > 4 does not 
matter, as the rules are the same. Scheme II cannot be considered as DRED or DREG, since 
the results obtained correspond to an anticommuting-75 scheme of DREG when there are an 
even number of 75's, but DREG gives no way to define traces with an odd number of 75's 
and maintain both cyclicity and anticommutativity of the 75. 



6 Summary 

To summarize, we have explicitly demonstrated that DREG leads, in a non-trivial fashion, 
to the correct factorized form for the transverse Drell-Yan squared amplitude in the limit of 
collinear gluon radiation. The resulting n-dimensional one-loop transversity splitting function 
is in agreement with that obtained in [2] as an intermediate step in the determination of the 
corresponding two-loop splitting function. The result is the same in the anticommuting-75 
scheme and in the HVBM scheme. This implies that the one-loop transversity distributions 
and their two-loop evolutions are the same in both schemes. The transversity distribution 
which would be obtained from the transverse Drell-Yan process using DREG subprocess 
cross sections was shown to be consistent with the corresponding one of DRED. In light of 
the comparative simplicity of the latter scheme, a minor variant on DRED was considered 
(scheme III) which can help in avoiding the addition of the DRED counterterms, at the very 
least. The link between DRED and the anticommuting-75 scheme of DREG was clarified and 
a specific anticommuting-75 scheme (scheme II) was formulated making use of it. 
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Fig. 1. Leading order contribution to transverse 
Drell-Yan. 




Fig. 2. Bremsstrahlung contributions to trans- 
verse Drell-Yan. 
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